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Abstract. For given integers a, b and j > 1 we determine the set 
R^^\ of integers n for which a" — 6" is divisible by nK For j = 1,2, 
this set is usually infinite; we determine explicitly the exceptional 
cases for which a, h the set R^f\^ (j = 1, 2) is finite. For j = 2, we use 

Zsigmondy's Theorem for this. For j > 3 and gcd(a,6) = 1, i?^"*! 
is probably always finite; this seems difficult to prove, however. 

We also show that determination of the set of integers n for 
which a" + 6" is divisible by can be reduced to that of -R), j,- 



1. Introduction 

Let a, h and j be fixed integers, with j > 1. The aim of this paper 
is to find the set R^^\ of all positive integers n such that divides 
aP' — IP. For j = 1,2,..., these sets are clearly nested, with common 
intersection {!}. Our first results (Theorems [1] and [2]) describe this set 
in the case that gcd(a, 6) = 1. In Section H] we describe (Theorem [T^ 
the set in the general situation where gcd(a, b) is unrestricted. 

Theorem 1. Suppose that gcd(a, 6) = 1. Then the elements of the 
set R^^l consist of those integers n whose prime factorization can be 
written in the form 

n = p\^p\^ ...p^; {pi<P2<---<Pr, allh>l), (1) 

where pi \ a"' — 6"' (i = 1, . . . , r), with rii = 1 and rii = p\^P2^ . . ■p'^'Si 
(z = 2,...,r). 

In this theorem, the ki are arbitrary positive integers. The result is 
essentially contained in [TU], which described the indices n for which 
the generalised Fibonacci numbers Un are divisible by n. However, we 
present a self-contained proof in this paper. 

On the other hand, for j > 2, the exponents ki are more restricted. 

Theorem 2. Suppose that gcd{a,b) = 1, and j > 2. Then the elements 
of the set R^Jl consist of those integers n whose prime factorization can 
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be written in the form (QP, where 



Pi ^^''^ divides 




if Pi > 2; 
-b,a + b) ifpi = 2, 

and p?"^^^" I a"- - 6"- , with rii = p\^pl^ . . . p^^i (i = 2, . . . , r) . 

Thus we see that construction of n G -R^"'^ depends upon finding a 

prime pt not used previously with oT'^ —6"' being divisible by p{~^ . This 

presents no problem for j = 2, so that R^^l^ ^ ^^^^ ^^a\i ^'^^ usually 
infinite. See Section [5] for details, including the exceptional cases when 
they are finite. However, for j > 3 the condition p':~^ \ a"^' — 6"' is 

only rarely satisfied. This suggests strongly that in this case R^^l is 
always finite for gcd(a, b) = 1. This seems very difficult to prove, even 
assuming the ABC Conjecture. A result of Ribenboim and Walsh |9j 
implies that, under ABC, the powerful part of a" — 6" cannot often 
be large. But this is not strong enough for what is needed here. On 
the other hand, R'fl (j > 3) can be made arbitrarily large by choosing 
a and b such that a — 6 is a powerful number. For instance, choosing 
a = l+{qiq2 ■ ■ ■ QsY'^ and b = 1, where qi,q2, ■ ■ ■ ,qs are distinct primes, 
then R^^l contains the 2^ numbers (fx ^2 • • • " where the are or 1. 
See Example [7] in Section [71 

In the next section we give preliminary results need for the proof 
of the theorem. We prove it in Section El In Section jH we describe 
(Theorem [15]) i?^"'], where gcd(a, 6) is unrestricted. In Section [5l we 

find all a, 6 for which i?^^^ is finite (Theorem [TBI) . In Section [HI we 
discuss the divisibility of a" + 6" by powers of n. In Section [71 we give 
some examples, and make some final remarks in Section [HI 

2. Preliminary results 

We first prove a version of Fermat's Little Theorem that gives a little 
bit more information in the case x = \ (mod p\ 

Lemma 3. For a; G Z and p an odd prime we have 

^p-i + a;P-2 + ... + ^ + l^ (modp2) = l (modp); 

I 1 (mod p) otherwise . 

Proof. If X = 1 (mod p), say x = 1 + kp, then x^ = 1 + jkp (mod p^), 
so that 

p-i 

x^~^ + x^'"^ + -- - + x + l= p + kp''^^j = p (mod p^). (3) 

j=0 

Otherwise 

x(x-l)(x^~^H hx + l) = x^-x = (modp), (4) 
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SO that for x ^ 1 (mod p) we have + • • ■ + x + 1) = (mod p), 

and hence 

x^"^ + a;^"^H hx + l = x(x^"^H hx + l) + l = l (mod p). (5) 

□ 

The following is a result of Birkoff and Vandiver [2j, Theorem III]. 
It is also special case of Lucas [HI P- 210], as corrected for p = 2 by 
Carmichael [21 Theorem X]. 

Lemma 4. Let gcd(a, b) = 1 and p be prime with p \ a — b. Define 
t > by p^\\a — b for p > 2 and 2*|| lcm(a — b,a + b) if p = 2. Then for 
i>0 

pW\\aP'-bP\ (6) 
On the other hand, if p \ a — b then for £ > 

p\aP'-bP\ (7) 

Proof. Put X = a/b. First suppose that p is odd and p^\\a — b for some 
t > 0. Then as gcd(a, b) = 1, b is not divisible by p, and we have x = 1 
(mod p*). Then from 

aP-bP = (a- b)bP-\xP-' + x^-^ + ■ ■ ■ + X + 1) (8) 

we have by Lemma [S] that p*+^||a^ — bP. Applying this result £ times, 
we obtain ([6]). 

For p = 2, we have p*"'"^||a^ — b^ and from = 6^ = 1 (mod 8), we 
obtain 2^||a^ + 6^, and so p^+^Ha^ — 6^. An easy induction then gives 
the required result. 

Now suppose that p \ a — b. Since gcd(a, &) = 1, (I7j) clearly holds if 
p I a or p I 6, as must happen for p = 2. So we can assume that p is 
odd and p \ b. Then x ^ 1 (mod p) so that, by Lemma [3] and ([H]), we 
have p\ aP — bP. Applying this argument £ times, we obtain ([7]). □ 

For n G R^l, we now define the set V^^l{n) to be the set of all prime 

powers p'^ for which np^ G R^P^- next result describes this set 

precisely. (Compare with [10, Theorem 1(a)]). 

Proposition 5. Suppose that j > 1, gcd(a, 6) = 1, n G R^^l and 

_2'=2 JJ/^', n = JJp^'' (9) 

p>2 p 

and define 62 by 2*^^ || lcm(a'^ — b^, a" + 6") . Then 

p(i)(n)= [j {pf'^keN}, (10) 

p|a"— 6" 

and for j > 2 

r!f,lin)= U |/:l<A;<L^^^j|. (11) 
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(2) 

Note that 62 is never 1. Consequently, if 2m & R\l, where m is odd, 
then 4m G R^al- Also, 2 G R^^l for j < 3 when a — 6 is even. 

Proof. Taking n G R^^\ we have, from ([9]) and the definition of 62 that 
jkp < Cp for all primes p. Hence, applying Lemma H] with a, h replaced 
by a" , 6" we have for p dividing aJ^ — 6" that for £ > 

pep+£||^n/ (12) 

So (np^)-^ I a"''^ — is equivalent to i{kp + i) < ep + £, or (j — l)i < 
— jkp. Thus we obtain (ITOl) for j > 2, with £ unrestricted for j = 1, 
giving (fTUjl . 

On the other hand, if j9 j" a" — 6", then by Lemma H] again, j 
_ 5V ^ so that certainly {npy \ a'^p' - □ 

We now recall some facts about the order function ord. For m an 
integer greater than 1 and x an integer prime to m, we define ordm,(a;), 
the order of x modulo m, to be the least positive integer h such that 
x^ = 1 (mod m). The next three lemmas, containing standard mate- 
rial on the ord function, are included for completeness. 

Lemma 6. For x G N and prime to m we have m \ x"' — 1 if and only 
if oidmix) I n. 

Proof. Let ordm(a;) = h, and assume that m \ x" — 1. Then as m \ x^—1, 
also m I x^^'^^^'"'^ — 1. By the minimality of h, gcd(/i, n) = h, i.e., h \ n. 
Conversely, ii h\n then x'^ — 1 | x"' — 1, so that m | x"' — 1. □ 

Corollary 7. Let j > 1. We haven^ \ x"— 1 if and only ifgcd{x,n) = 1 
and ord„j (x) \ n. 

Lemma 8. For m = YlpP^'' and x G N and prime to m we have 

ordm(x) = lcmpordpfep(x). (13) 

Proof. Put hp = ordp/p(x), h = ordm(a;) and h' = Icmphp. Then by 
Lemma [6] we have p-^^ | x'*' — 1 for all p, and hence m | x^' — 1 . Hence 
h I h'. On the other hand, as p-^^ \ n and m | x'^ — 1, we have p-^^ \ x^ — 1, 
and so hp \ h, by Lemma El Hence h' = Icm^ hp \ h. □ 

Now put = ordp(x), and define t > by ]9*||x^* — 1. 

Lemma 9. For gcd(x, n) = 1 and i > we have | p — 1 and 
ordp.(x) =p™^-(^-*'0)p,. 

Proof. Since p \ x^~^ — 1, we have p* | p — 1, by Lemma El Also, from 
p^ I x'^^'^p"^^^ — 1 we have p | x^^'^p'^^^ — 1, and so, by Lemma El again, 
p^ = ordp(x) I ordp«(x). Further, if i < t then from p^ | x^* — 1 we have 
by Lemma El that ordp«(x) | p*, so ordp«(x) = p^. Further, by Lemma 
H for It > t 

_ 1^ (14) 
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SO that, taking u = i > t and using Lemma El oidpt^x) \ 

Also, if t < u < i, then, from (fT^ . x^* ^ 1 (modp^). Hence 

ordp^" (x) = p^~^p^ for £ > t. □ 

Corollary 10. Lei j > 1. For n = YlpP^'' '^'^^ x G N and prime to n 
we have \ — 1 if and only if gcd{x, n) = 1 and 

icmp/pp* I n^''"- ^^^^ 

p 

Here the kp = max{jkp — tp, 0) are integers with tp > 0. 

Note that p*, kp and tp in general depend on x and j as well as on p. 
What we actually need in our situation is the following variant of 
Corollary [ini 

Corollary 11. Let j > 1. For n = YlpP'''' and integers a,b with 
gcd(a, b) = 1 we have \ a^—h^ if and only ifgcd{n, a) = gcd(n, b) = 1 
and 

icuipp^pp^ I n^'""- ^^^^ 

p 

Here the kp = max(jfcp — tp, 0) are integers with tp > 0. 

This corollary is easily deduced from the previous one by choosing x 
with bx = a (mod n^). 

By contrast with Proposition El our next proposition allows us to 
divide an element n G R^^l by a prime, and remain within R^l- 

Proposition 12. Let n G R^^\ with n > 1, and suppose that Pmax is 
the largest prime factor of n. Then n/pmax G R^^l- 

Proof. Suppose n G R^^b^ so that (fT^ holds, with x = a/b, and put 
q = Pmax- Then, since for every p all prime factors of are less than 
p, the only possible term on the left-hand side that divides q'^'' on the 
right-hand side is the term q^i. Now reducing kq by 1 will reduce k'^ 
by at least 1, unless it is already 0, when it does not change. In either 
case ( fT5|) will still hold with n replaced by n/q, and so n/g G R^^l- □ 

Various versions and special cases of Proposition [T2] for j = 1 have 
been known for some time, in the more general setting of Lucas se- 
quences, due to Somer [HI Theorem 5(iv)], Jarden P Theorem E], 
Hoggatt and Bergum [5], Walsh [13], Andre- Jeannin [1] and others. 
See also Smyth [ini Theorem 3]. 

In order to work out for which a, b the set R^^l is finite, we need the 
following classical result. Recall that a" — 6"" is said to have a primitive 
prime divisor p if the prime p divides a"^ — b"^ but does not divide — b^ 
for any k with 1 < A; < n. 
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Theorem 13 (Zsigmondy [TT]). Suppose that a and b are nonzero 
coprime integers with a > b and a + b > 0. Then, except when 

• n = 2 and a + b is a power of 2 
or 

• n = 3, a = 2, b = —1 

or 

• n = 6, a = 2, b = 1, 

a" — b^ has a primitive prime divisor. 

(Note that in this statement we have allowed b to be negative, as did 
Zsigmony. His theorem is nowadays often quoted with the restriction 
a > 6 > and so has the second exceptional case omitted.) 

3. Proof of Theorems [T] and [2] 

Let n G R^^l have a factorisation ([T]), where Pi < P2 < ■ ■ ■ < Pr and 

all ki > 0. First take j > 1. Then by Proposition [T2] n/ p^'' = G R^l, 
and hence 

("-/PrO/Pr-? = ^r-l, ■ ■ ■ , Pi' = Ui, 1 = Ui 

are all in R^l- Now separate the two cases j = 1 and j > 2 for 
Theorems [T] and [2] respectively. Now for j = 1 Proposition [5] gives us 
that Pi I a"' — b""^ [i = 1, . . . , r), while for j > 2 we have, again from 
Proposition [5], that 



a — b if pi > 2; 

lcm(a — b,a + b) if pi = 2, 



pU 1)^1 (^livides 

and I a"' — b""^ {i = 2, . . . , r). Here we have used the fact that 

gcd{pi,ni) = 1, so that if p^' \ (a"" - b''*)/n'f then p'l' \ a"* - b""' (i.e., 
we are applying Proposition [5] with all the exponents kp equal to 0.) 

4. Finding i?^-'] when gcd(a, 6) > 1. 

For a > 1, define the set JF^ to be the set of all n G A/" whose prime 
factors all divide a. To find R^^l in general, we first consider the case 
b = 0. 

Proposition 14. We have -R^^q = -^i^o — -^a, while for j > 3 the set 

^al = Sa \ where Sa '' is a finite set. 

Proof. From the condition \ a"^, all prime factors of n divide a, so 
rI^I C J^a, say rI^I = J^a\ Si^\ We need to prove that Si^^ is finite. 
Suppose that a = pi' . . .p"'', with pi the smallest prime factor of a. 
Then n = p^' . . . p^"^ for some fc, > 0. From \ a" we have 

A:.<-P^...p^ (^ = l,...,r). (17) 

J 
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For these r conditions to be satisfied it is sufficient that 

E*^ , mm^_^ a,- t'' , , 

k< t^pp-'^'. 18 

Now (|T8l) holds if j = 1 or 2, as in this case, from the simple inequality 
k < 2^^^^ valid for all G N, we have 

j2h<h^U^ <^^^^^p^--^'\ (19) 
i=i ^ 

Hence Sa'^ is empty if j = 1 or 2. 

Now take j > 3, and let K = Ka^ be the smallest integer such that 
Kp:[^ < (min[^ia,)/j. Then holds for EI=i > and S'^^^ is 
contained in the finite set S" = {n G N, n = . . p^'' : X]i=i < 
(To compute S'i"''' precisely, one need just check for which r-tuples 
(/ci, . . . , kr) with ki < K any of the r inequalities of (|T7j) is vio- 
lated. □ 

One (at first sight) curious consequence of the equality R^^\ = R^^\ 
above is that n \ a" implies | a". 

Now let g = gcd(a, b) and a = aig, b = big. Write n = Grii, where 
all prime factors of G divide g and gcd{ni,g) = 1. Then we have the 
following general result. 

Theorem 15. The set R^^l is given by 
Ral = {n = Gm:Ge Tg, G ,G and gcd((7, m) = 1} \ i?, (20) 

•where R is a finite set. Specifically, alln = Gni G R have 1 < ni < j /2 
and 

G = qi\..qt, (21) 

where 

m 

Y,^.<K%. (22) 



i=l 



Here the qi are the primes dividing g, and K^JX is the constant in the 
proof of Proposition [I^] above. 

Proof, supposing that n G i?^^^ we have 

I a" - 6" (23) 
and so \ g^{a^ — 6"). Writing n = Gni, as above, we have 

n{ I (af )"i - (6f )"i (24) 

and 

G^ I ((af )"' - (&f )"') • (25) 
Thus ( l23l) holds with n, a, b replaced by ni, af, bf. So we have reduced 
the problem of (!23|) to where gcd(a, b) = 1, which we can solve 
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for rii prime to g, afong with the extra condition Now, from the 
fact that R^gl = Tg from Proposition [TH we have \ and hence 
I fif"^"! for all G e , provided that ni > j/2. Hence ([25]) can fail 
to hold for all G E J-'g only for 1 < rii < j/2. 

Now fix ni with 1 < rii < j/2. Then note that by Proposition [HI 
Qj I gGni g^j^j hence fl23|) holds for all G G \ S, where S* is a finite 
set of G's contained in the set of all G"s given by fl^ and ([22]) • □ 

Note that (taking ni = 1 and using (1^ ) we always have R^Jl C i?^"'^. 
See example [7] in Section [3 

5. When are R^^j^ and _R^^^ finite? 

First consider -R^^^. From Theorem [T] it is immediate that R^^l con- 
tains all powers of any primes dividing a — b. Thus R^^l is infinite unless 

a — 6 = lb 1, in which case R^^l = {!}■ This was pointed out earlier by 
Andre- Jeannin [H Corollary 4]. 

(2) 

Next, take j = 2. Let us denote by I the set of primes that divide 
some n E R\ l and, as before, put g = gcd(a, b). 

Theorem 16. The set R^^\ = {1} if and only if a and b are consecutive 

(2) (2) 

integers, and Rah = {1,3} if and only if ab = —2. Otherwise, R^l is 
infinite. 

If R^^^^^l^ = {1} (respectively, = then V^^l '^^ ^/ 

prime divisors of g (respectively, 3g). Otherwise V^^l ^■^ infinite. 

The application of Zsigmondy's Theorem that we require is the fol- 
lowing. 

Proposition 17. If R\ l contains some integer n > 4 then both R\ ^^ 
and V^al infinite sets. 

Proof. First note that if a = 2, 6 = 1 (or more generally a — b = ±1) 
then by Theorem [21 R^"^^ = {1}. Hence, taking n G R^^l with n > 4 
we have, by Zsigmondy's Theorem, that a" — 6" has a primitive prime 
divisor, p say. Now if p \ n then, by applying Proposition [12] as many 

(2) 

times as necessary we find p \ n', where n' G R^l, and now p is the 
maximal prime divisor of n'. Hence, by Proposition [T2] again, n" = 
n'/p G R^al from n' = pn" and Proposition [5] we have that 

p I a"'" — h^" , contradicting the primitivity of p. 

Now using Proposition [5] again, np G R^al- Repeating the argument 
with n replaced by np and continuing in this way we obtain an infinite 
sequence n, np, nppi, nppiP2, . . . , nppiP2 . . .pi, . . . of elements of R^^h 
where p < pi < p2 < ■ ■ ■ < Pi < ■ ■ ■ are primes. □ 
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Proof of TheoremlT^ Assume gcd(a, 6) = 1, and, without loss of gen- 
erality, that a > and a > b. (We can ensure this by interchanging a 
and b and/or changing both their signs.) If a — 6 is even, then a and b 
are odd, and — 6^ = 1 (mod 2*"''^), where t > 2. Hence 4 G R^^h 

Proposition El and so both R^^l and V^l are infinite sets, by Proposi- 
tion [HI 

If a — 6 = 1 then R^"^^ = {!}, as we have just seen, above. 
If a — 6 is odd and at least 5, then a — b must either be divisible by 
9 or by a prime p > 5. Hence 9 or p belong to R\l, by Proposition [5], 

and again both R^^j^ and V^l are infinite sets, by Proposition [TTl 

If a - 6 = 3 then 3 G i?^, and - b^ = 9(6^ + 3b + 3). If 6 = -1 
(and a = 2, ab = —2) or —2 (and a = 1, ab = —2) then — b^ = 9 and 
so, by Theorem [21 so i?^^^ = {1,3}. Otherwise, using gcd(a, 6) = 1 we 
see that — b^ > 5, and so the argument for a — b > 5 but with a, b 
replaced by a^, b^ applies. □ 

6. The powers of n dividing a" + 6" 

Define i?^-^]^ to be the set {n G N : divides a"" + 6"}. Take j > 1, 
and assume that gcd(a, 6) = 1. (The general case gcd(a, 6) > 1 can be 
handled as in Section HI) We then have the following result. 

Theorem 18. Suppose that j > 1, gcd(a, 6) = 1, a > and a > \b\. 
Then 

(a) R^^l^ consists of the odd elements of R^^\^^, along with the num- 
bers of the form 2ni, where ni is an odd element of R^^^ 

(b) If j > 2 the set R^l^ consists of odd elements of R^f^_^ only . 

Furthermore, for j = 1 and 2, the set R^^l^ is infinite, except in the 
following cases: 

• If a + b is 1 or a power of 2, (j, a, b) ^ (1, 1, 1), when it is {!}; 
. <J+ = {1,2}; 

. 4'i+ = {l,3}. 

Proof. If n is even and j > 2, or if 4 | n and j = 1, then | a" + 6" 
implies that 4 | a" + 6", contradicting the fact that, as a and b are not 
both even, a" + 6" = 1 or 2 (mod 8). So either 

• n is odd, in which case \ a" + 6" is equivalent to finding the 

odd elements of the set R^a\; 
or 

• j = 1 and n = 2ni, where ni is odd, and belongs to R^^2 _i,2- 
Now suppose that j = 1 or 2. If a + 6 is ±1 or ± a power of 2, then, 
by Theorem [2l all n G R^P_i, with n > 1 are even, so for j = 2 there are 
no n > 1 with | a" + 6" in this case. Otherwise, a + b will have an 
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odd prime factor, and so at least one odd element > 1. By Theorem [TBI 

(2) 

and its proof, we see that R^^ ij, will have infinitely many odd elements 
unless a{—b) = —2, i.e. a = 2, 6 = 1 (using a > and a > \b\). 

For j = 1, there will be infinitely many n with n \ a"^ + b"', except 
when both a + b and + 6^ are 1 or a power of 2. It is an easy exercise 
to check that, this can happen only for a = 6 = 1 or a = 1, 6 = 0. □ 

If g = gcd(a, b) > 1, then, since R^^l^ contains the set -R^'^O) "^i^^ ^e 
infinite, by Proposition [TH For j > 3 and gcd(a, 6) = 1, the finiteness 
of the set R^l^ would follow from the finiteness of R^l, using Theorem 

mh)- 

7. Examples. 

The set R^^l has a natural labelled, directed-graph structure, as fol- 
lows: take the vertices to be the elements of R^Jl, and join a vertex n to 

a vertex np as n -^p np, where p G Vj^^l . We reduce this to a spanning 
tree of this graph by taking only those edges n —>-p np for which p is 
the largest prime factor of np. For our first example we draw this tree 
(Figure [T]). 

1. Consider the set 

Rf} =1, 2, 4, 20, 220, 1220, 2420, 5060, 13420, 14740, 23620, 55660, 
145420, 147620, 162140, 237820, 259820, 290620, 308660, 
339020, 447740, 847220, 899140, 1210220, . . . , 

(sequence A127103 in Neil Sloane's Integer Sequences website). 
Now 

- 1 = 2^ ■ 5^ ■ 11^ • 61 • 1181, 
showing that ^£^20) = {11, 11^61, 1181}. Also 
3220 _ 1 = 2^ ■ 5^ ■ 11^ • 23 ■ 61 ■ 67 ■ 661 ■ 1181 ■ 1321 ■ 3851 ■ 5501 
■ 177101 ■ 570461 ■ 659671 ■ 24472341743191 ■ 560088668384411 
■927319729649066047885192700193701, 

so that the elements of ^3^^(220) less than 10V220, needed for 
Figure [H are 

11,23,61,661,1181,1321,3851. 

2. Now 

= 1, 2, 3, 4, 6, 12, 21, 42, 52, 84, 156, 186, 372, . . . , 
whose odd elements give 

^ 3^ 21, 609, 903, 2667, 9429, 26187, .... 
See Section [61 
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55660 147620 162140 237820 30S660 339020 



Figure 1 . Part of the tree for i?3 {, showing all elements 
below 10^ 

3. We have 

4S^ = <-2 = 1,5,55,1971145,..., 

(2) 

as all elements of _2 are odd. Although this set is infinite by 
Theorem [TU| the next term is 1971145p where p is the smallest 
prime factor of 3^^^^^^^ + 2^^^^^^^ not dividing 1971145. This 
looks difficult to compute, as it could be very large. 

4. We have 

<-3 = = 1,7,2653,.... 

Again, this set is infinite, but here only the three terms given 
are readily computable. The next term is 2653p where p is the 
smallest prime factor of 4^^^^ + 3^^^^ not dividing 2653. 

5. This is an example of a set where more than one odd prime 
occurs as a squared factor in elements of the set, in this case 
the primes 3 and 7. Every element greater than 9 is of one of 
the forms 21m, 63m, 147m, or 441m, where m is prime to 21. 

Rfl^ =1, 3, 9, 21, 63, 147, 441, 609, 1827, 4137, 4263, 7959, 

8001, 12411, 12789, 23877, 28959, 35931, 55713, 56007, 
86877, 107793, 119973, 167139, 212541, 216237, 230811, 
232029, 251517, 359919, 389403, . . . , 

6. i^^yjjoi,! = {1, 2, 3, 5, 6, 10, 15, 30}. This is because 27001 - 1 = 
2^ ■ 3^ ■ 5^, and none of 2700 1*^ — 1 has a factor p^ for any prime 
p > 5 for any n = 1, 2, 3, 5, 6, 10, 15, 30. 

7. Rfl^ = {1,2,3,6,42,1806}? Is this the entire set? Yes, unless 
X9I8O6 _ ]^ ig divisible by p"^ for some prime p prime to 1806, 
in which case 1806p would also be in the set. But determining 
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whether or not this is the case seems to be a hard computational 
problem. 

8. R'^q 2^ example with gcd(a, 6) > 1. It seems highly probable 
that 

2 = (^2\{2,4,8})U(3^2) 
= 1, 3, 6, 12, 16, 24, 32, 48, 64, 96, 128, 192, 256, 384, 512, 768, 1024, . 

However, in order to prove this. Theorem [13 tells us that we 
need to know that 28^ ^ 1 (mod p^) for every prime p > 3 
and every £ > 0. This seems very difficult! Note that = 
J-2\{2,4,8} and i?gi = {l,3}. 



8. Final remarks. 

1. By finding R^l, one is essentially solving the exponential Dio- 
phantine equation x^y = — b^, since any solutions with x < 
are readily found. 

2. It is known that 

R^^h = |n G N : n divides — |. 

a — b 

See [TOl Proposition 12] (and also Andre- Jeannin [1], Theorem 
2] for some special cases of this result.) This result shows that 
r[^1 = {n G N : n divides where the Un are the generalised 
Fibonacci numbers of the first kind defined by the recurrence 
uq = 1, ui = 1, and Un+2 = (a + b)un+i — abun {n > 0). This 
provides a link between Theorem [1] of the present paper and 
the results of [TUj . 

The set R^^]^ is a special case of a set {n G N : n divides t;„}, 
also studied in [lO]. Here {vn) is the sequence of generalised 
Fibonacci numbers of the second kind. For earlier work on this 
topic see Somer [T^ . 

3. Earlier and related work. The study of factors of a" — 6" dates 
back at least to Euler, who proved that all primitive prime fac- 
tors of a" — fe"" were = 1 (mod n). See j2| Theorem 1]. Chapter 
16 of Dickson ^ (Vol 1) is devoted to the hterature on factors 
of a"" ± 6". 

More specifically, Kennedy and Cooper [7] studied the set 
R^iQi- Andre- Jeannin fTj Corollary 4] claimed (erroneously - 
see Theorem [T8|) that the congruence + b^ = (mod n) 
always has infinitely many solutions n for gcd(a, b) = 1. 

4. Acknowledgement. I thank Hugh Montgomery for telling me 
about Zsigmondy's Theorem. 
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